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Abstract. Methods for estimation of different reference times which appear in the description of 

transition of a strong adiabatic shock into the radiative era are reviewed. The need for consideration of an 

additional transition subphase in between the end of the adiabatic era and the beginning of the radiative 

fvj . "pressure-driven snowplow" stage for a shock running in the uniform or nonuniform medium is emphasized. 

, , ' This could be of importance in particular for studying of the interaction of supernova remnants (SNRs) 

^ |. with molecular clouds and therefore for understanding the processes of the cosmic ray production in such 

■^C ' systems. The duration of this subphase - about 70% of SNR age at its beginning - is almost independent 

of the density gradient for media with increasing density and is longer for higher supernova explosion 

energy and for smaller density in the place of explosion. It is shown as well that if the density of the 

ambient medium decreases then the cooling processes could differ from the commonly accepted scenario 

^ ' of the "thin dense radiative shell" formation. This property should be studied in the future because it is 

C""~~ , important for models of nonspherical SNRs which could be only partially radiative. 
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Physical processes accompanying the evolution of supernova remnants (SNRs) is a complex system. 
It is almost impossible to account for all of them in a single model of SNR. Therefore, the whole 

O I evolution of SNR from a supernova explosion until the mixing of a very old object with the 

interstellar matter is divided on a number of the model phases (e.g. 0121121): the free-expansion, 

^ I adiabatic, radiative and dissipation stages. There are some physical processes important during a 

given stage, some others could be neglected. Such an approach allows for rather simple analytical 
description of SNR evolution during each phase. 
S^ i The role of radiative losses, which is negligible in the adiabatic phase of SNR evolution, becomes 

H I more and more promiment with time. They are so important in old SNRs, that they essentially 

modify the dynamics of such SNRs. Theoretical systematization of timescales and the role of 
different physical process in cooling of adiabatic SNR was first reviewed in 0]. In general, the 
transition to the radiative stage can be studied numerically, by following the history of the shocked 
flow as it is done e.g. by [S1|S1|71|S1 El El ■ The analytical treatments are of great importance as 

well, e.g. [niiiiiniiiiiisiiini- 

The physical processes in the radiative blast wave, namely, quick cooling of an incoming flow 
and formation of the thin dense cold shell which moves due to the pressure of internal gas makes 
the so called "pressure-driven snowplow" (PDS) model within the "thin-layer" approximation to 
be adequate for description of this stage of SNR evolution [TTl [T51 ITSl . 

The PDS model was introduced by McKee & Ostriker; their analytical solution JEIEI widely 
used for the description of evolution of the radiative shell gives a power-law dependence i? ex i™ 
(where t is age and i? is a position of the shock) with constant m (which equals to 2/7 for the 
uniform medium). However, numerical studies cited above give something different values of the 
deceleration parameter m (defined as to = d\nR/dlnt), namely ~ 0.33 (51^1. We have shown 
analytically in ^^1 that the evolution of the radiative shell is given by variable m and that the 
discrepancy between the analytical and numerical results is only apparent. In fact, the usage 
of McKee & Ostriker analytical solution assumes that SNR has already reached the asymptotic 
power-law regime with constant value oi m — 2/7. The time needed to reach this asymptotic 
regime is however long compared to the SNR age. 

It is common for an approximate theoretical description of SNR evolution to simply switch 
from the adiabatic solution to the PDS radiative one at some moment of time. However, we stress 



in this paper the result visible also in previous calculations, namely, the need for an intermediate 
transition subphase between the adiabatic and radiative stages, with duration more than a half 
of SNR age it has at the time when radiative losses of gas passing through the shock begins to 
be prominent. Thus the radiative era which begins after the end of the adiabatic one, have to 
be divided on two phases: the transition subphase, when the radiative losses become to modify 
dynamics and to lead to the formation of the thin radiative shell, and the PDS stage when one can 
apply the PDS analytic solution. In the present paper the role of nonuniform interstellar medium 
on the duration of the transition subphase is considered. 

2 Transition to the radiative phase 

2.1 Definitions of different reference times 

Let us consider the spherical shock motion in the medium with the power-law density variation 
Po{R) = AR~'^ , where A and cj are constant; indexes "o" and "s" refer hereafter to the pre- 
and post-shock values. The dynamics of the adiabatic shock in such a medium is given by Sedov 
solutions 1^ where the shock velocity D oc i?-(3-'^)/2 and R oc i2/(5-t^)_ 

Moving through medium, the shock decelerates if the ambient density distribution increases or 
does not quickly decrease {uo < 3). The shock temperature Tg oc D^ decreases with time as well. 
Starting from some age tiow when Tg = Tjow ~ 3 x 10^ K, which corresponds to the minimum of 
the cooling function A(T), the radiative losses of shocked plasma are more and more prominent 
with falling of T (Fig. ^. The maximum in the energy losses is when the shock temperature 
Tg = Thi ~ 2 X 10^ K, the corresponding Sedov time (i.e. calculated under the assumption that the 
shock is adiabatic up to this time) is thi. 

There is a number of reference times in between tiow and thi ^ 1201 IIT)] . Once a parcel of gas is 
shocked its temperature changes due to expansion and cooling Ta = Ta,oxp + ^l^^rad, where the dot 
marks the time derivative. One may define the "dynamics-affected" time tdyn by the equation 

-^ a,cxpv''dyn/ -^ a,rad\^dynj ■ \^} 

If a fluid element is shocked after this time, its temperature decreases faster due to radiation 
than as a consequence of expansion. At other time tsag, the radiative cooling begins to affect the 
temperature distribution inside the shock. When the rate of change of the shock temperature Tg 
begin to be less than Tai the temperature downstream of the shock will sag rather than rise. Thus 
the equation for tgag is 

-'sl^t^sagj — -'at.t^sagj. \'^) 

Radiative losses cause the faster - comparing to the adiabatic phase - deceleration of the 
forward shock. This faster deceleration begins to be prominent around the "transition age" ttr 
when the shock pressure decrease due to the radiative losses becomes to be effective. Then, the 
shocked gas radiates away its energy rather quickly, cools till the temperature T '^ 10"* K and 
forms a dense shell. The formation of the shell is completed around the "time of shell formation" 
tgf which is larger than ttr! the latter which marks the end of adiabatic era. After tgf the thermal 
energy of all swept-up gas is rapidly radiated and the thin dense shell expansion is caused by the 
thermal pressure of the interior. 

The time tiow is given by the equation 

TsiUovj) = Tiow (3) 

A similar equation defines the time thi 

Tg(thi) = Thi, (4) 

which was suggested to be a measure of ttr ^lO]' However, as we shall demonstrate later, the 
post-shock temperature of plasma at ttr is of order 10^ K > Thi and thi is larger than ttr in about 
3.5 times (Sect.|OJ- Therefore it is not correct to calculate the "highest-losses" SNR age with the 
shock motion law valid during the adiabatic era. 

A simple approach to locate ttr bases on the comparison of the radiative losses with the initial 
thermal energy of the shocked fluid JOI- A shocked fluid element cools during the cooling time 



Atcooi oc e(Ts, ps)/A{Ts, Ps), where e = (7 — 1)~^ psk-QTs/ purn^ is its initial thermal energy density, 7 
is the adiabatic index, k^ is the Boltzman constant, mp is the proton mass. During the adiabatic 
phase the cooling time is larger than SNR age t. The radiative losses may be expected to modify 
dynamics when the cooling time Aicooi < i- In such approach the transition time is a solution of 
equation 

itr = Aicool(itr)- (5) 

Let us assume that the cooling function A ex n^T^^ with /3 > and n is the hydrogen number 
density, then Aicooi ot n~^T^~^^ oc t~6(i+'')/^ with the use of Sedov solutions for uniform medium. 
For the shock running in the power-law density distribution, the upstream hydrogen number density 
and the post-shock temperature at time t is 

no oc i-2.^/(5-"), T, oc t-2(3-c.)/(5-^)_ ^g) 

Therefore Aicooi oc t^'^ with i] ~ (2(3 — uj){1 + P) — 2uj)/{5 — uj)) for such density distribution. For 
P — 1/2 the index 77 is the same as found in {QJ. 

The way to estimate the time of the shell formation tgi was suggested in [2011^. If an element 
of gas was shocked at time is then the age of SNR will be t^ = ts + Aicooi(is) when it cools down. 
The minimum of the function tdtg) has the meaning of SNR age when the first element of gas 
cools and is called "SNR cooling time" icooi- Let ii be the time when the shock encountered the 
fluid element which cools first. If so, tc — ti{ts/ti) + Atcooi{ti){ts/ti)^^ . Setting dtc/dts\t^=tx — 
one obtain that 

icool= (l+?7)Aicool(il), (7) 

icooi 1 + ?7 /o^ 

-7 — ; \^) 

The cooling time tcooi > ii by the definition, therefore it must be that 77 > 0. This is the case for 

c^<3(l + /3)/(2 + /3); (9) 

that is w < 2 (9/5) for /3 = 1 (1/2). The equation 

ti=?7Aicooi(ii) (10) 

is more suitable for practical use than {T)). If the medium is uniform then icooi = 17<i/12 for /3 ~ 1 
and icooi = 14ii/9 for /3 = 1/2. 

The "SNR cooling time" icooi = min(tc) was initially suggested to be taken as the time of the 
shell formation. Numerical experiments for shock in the uniform medium suggest that isf is a bit 
higher (of order 10%) than icooi |22] and the reason of this could be that the compression of the 
shell is also effective after cooling of the first element and takes additional time. 

Another point is that the solution for adiabatic shock used in © might not formally be ap- 
plicable there because ti > ttr (see Eq. ^^). We believe however that the level of accuracy in 
estimation of ttr, the small difference between ttr and ti (about 30% in the case of uniform medium, 
Sect. l3.l|l as well as close values of icooi and tgf allow one to use the Sedov solution in © and to 
assume tgi ~ tcooi- 

We would like to note once more that the transition time itr is an approximate estimation 
on the end of the adiabatic stage and beginning of the radiative era, while the time of the shell 
formation tsf marks the time when one can start to use the PDS model where hot gas pushes the 
cold dense shell^ . The structure of the flow re-structurises and the shell forms during the transition 
subphase given by the time interval (itr, tsf). We shall demonstrate later that the ratio isf/itr with 
itr given by (O and tgf by ||SJ) is always larger than unity (see Eq. IJSHI) and that the transition 
subphase is not short as it is generally assumed. 

One more time, namely the "intersection time" ti G (itrj^sf) was introduced in jIQ, as a time 
when two functions - the adiabatic dependence R = R{t) (valid before itr) and the PDS depen- 
dence Rsh — Rsh{t) (valid after tgf) - intersect being extrapolated into the transition subphase. 
This intersection time could be useful in some tasks when the level of accuracy is such that one 
may sharply switch from the adiabatic solution to the radiative one without consideration of the 
transition subphase. 

^The PDS analytical solutions which describe the evolution of SNR after the shell formation time are presented 
in II3III41II6I for uniform ISM and in 1141 for ISM with power-law density variation. 



2.2 Cooling time 

The expression 

_ £(rs,Ps) , . 

Aicool - -T-p^ 7 (11) 

used in jlOj to calculate the cooling time, equates the energy losses AAicooi with initial thermal 
energy density eg of a fluid element under condition that the density and temperature of this element 
are constant. More detailed model should account for the density and temperature history during 
Atcooi- Namely the above equation should be replaced with a differential one: 

dt/dt=-K{T,p). (12) 

The total internal energy U = eV oi gas within the volume V changes as dU = TdS — PdV where 

S is entropy and P is pressure. The evolution of the thermal energy per unit mass E — e/ p is 

therefore 

dE P fdp\ ds , , 

^ ^\ T- (13) 



dt p2 yQij Qi 

where s = (3fcB/2''np/i) In [P/ p"') is the entropy per unit mass {m-p is the mass of proton, p is the 
mean particle weight). So, Eq. \i2\ becomes 

t| ^ -^<M. (14, 

at p 

here the temperature T, density p, pressure P, energy E are functions of Lagrangian coordinate a 
and time t. 

As it follows from H14(l and the definition of s, the time Atcooi may be also defined as a time 
taken for the adiabat P/ p'' to fall to zero. Kahn [^H] have found an interesting result. Namely, if 

/3=^ (15) 

7-1 

(that is (3 ~ 1/2 for 7 = 5/3) then one can derive Atcooi from (|14(l independently of the density 
and temperature history: 

AfKahn _ <^{Ts, Ps) ,.„. 

(P + l)A(Ts,ps) 

It can be checked that the same solution may be obtained from H13|) - (|14ll for any /3 if one assume 
that the gas is not doing work during Atcooi that is equivalent to putting dp/dt = in ifT^ . 
However, the density of fluid is not expected to be constant. In such situation one should solve 
the full set of the hydrodynamic equations which can be performed only numerically, while we are 
interested in a rather simple analytical estimation on cooling time for a general f3. Therefore it is 
more suitable to use the estimation (|11|) for the cooling time which follows just from comparison 
of the radiative losses with the initial energy. We shall see later that such approach describe the 
shock dynamics rather well (Fig.|21). 

2.3 Equations for the reference times 

Let us write equations for ftr and tgt for shock in nonuniform medium. We assume hereafter (3 = 1. 
Note that all the rest formulae can easily be modified if one uses (3 which coincides with a value 
given by (|15|l ; namely, as it follows from comparison of (|16|l and lllll , T in (|17|l have to be simply 
divided by /3 + 1. 

If the cooling function for a fiuid is approximately A = CT^^UeUH, where C is a constant, 
then (dJ yields 

Aicooi = T^ where T = — ^2/^, (I7) 

no[R) Cp.['y + 1) 

Pe is the mean mass of particle per one electron in terms of the proton mass (i.e. p = p.^riem-p — 
linrrip). The transition time itr is a solution of equation |(SJ): 

t^r=T 'f' (18) 



where the dependencies Ts(i), R{t) are those vahd on the adiabatic phase. The time ii can be 
estimated from (|10|l : 

no(i?(ti)) ^ ^ 

Now the SNR cooUng time tcooi and the time of the shell formation igf « icooi is given by IJHI- 
The estimations for the transition and the shell formation times are somewhat different in the 
literature because of different ways used to find the cooling time Aicooi and to approximate the 
cooling function A(T). 

For the adiabatic shock the rate of change of the shock temperature is 

t . _?ii^5. ,20, 

5 — w t 

Close to the shock, the fluid temperature in Sedov solution |19| is approximately 

T{a) / a N -"('T.'^) 



where a is Lagrangian coordinate. The value of k is given by 

a dT{a) 



(S) "'. <-) 



T{a) da 



(22) 



a=R 



where T{a) is the profile from Sedov solutions. It is k = 1 — Suj/A for 7 — 5/3 (see Appendix). 
Now we may find that the temperature in a given fluid element a changes due to expansion as 

. ^ 2(3 - c. - /.) r(a) 

-'a.Gxp ^^ ^ , ■ V "^z 

■ 5 ~ Ul t 

The rate Ta,rad due to cooling follows from dE/dt = — A/p: 

Ta.rad = - ^^T-iyiH (a)T(a)-^. (24) 

7 + 1 

Now we have to compare the above rates at the time ig, i-e. at the time when the parcel of fluid 
was shocked. The coordinate a = R{ts) by the definition. Thus Eq. |^ rewrites: 

tdvn = Z Al-cooU'dynj- (^5) 

5 — w 
Similarly, the equation for ig^g follows from ^i 

2k 

— UJ 

As one can see, the most of reference times are given by the equations of the form 

U = i^Aicooi(i»), (27) 

where i* is a given reference time and K is corresponding constant. It may be shown that the 
solution of such equation may be found as 

U = /v i/^i+^^itr. (28) 

The Sedov radius of the shock at this time is i?, = i^2/((5-c^)(i+r,))^^^^ 

2.4 The cooling function 

There are two choices of (3 in the literature, namely 1 and 1/2. The first case is used for nonequi- 
librium cooling model |24j where the cooling function for plasma with solar abundance may be 
approximated as jl()| 

A = lO^^Vn-HT"^ erg cm^^ s~^ (29) 
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Figure 1: Equilibrium (line 1) \2t)\ and noncquilibrium (line 2) |24| cooling functions, used in the 
literature to study the transition of SNRs into the radiative phase, and approximations H3()|l (line 
3) and 1)29(1 (line 4). The equilibrium cooling function from [21] is also shown for comparison (line 
5). 

This approximation is valid for range of temperatures T = (0.2 — 5) x 10^ K which is important 
for description of transition into the radiative phase. Another possibility is to use the equilibrium 
cooling model as it was done in [2S1 1201 M I2H1 122] • In this case the approximate proportionality 
A oc T~^/^ is a reasonable one, e.g. for results on the cooling of the collisional equilibrium plasma 
from [2nilS]j the actual approximation 



A = 1.3 X 10~'^nenffT 



^1/2 



erg cm s 



(30) 



is written for plasma with almost the same abundance as above and is valid for T = (0.05 — 50) x 

lo^K i2n]. 

Different cooling functions are compared with their approximations on Fig. ^ At lower tem- 
peratures, the nonequilibrium cooling is less effective in energy losses than the equilibrium one 
(compare lines 2 and 5). This is because the cooling rate for temperatures higher than ~ 3 x 10^ K 
is mostly due to free-free emission while below this temperature the cooling is mostly due to the 
line emission from heavy elements (most heavy elements are completely ionized above ~ 3 x 10^ K). 
Under nonequilibrium ionization conditions the ions are underionized because electrons are much 
colder than ions and thus there is less emission from ions [S1|2H] (see also Fig. 18 in 24 ). 



3 Reference times and transition subphase 

3.1 Shock in a uniform ISM 

Let us compare the sequence of different reference times with numerical calculations 10 of transi- 
tion of the adiabatic shock into the radiative era, on example of the shock motion in the uniform 
ambient medium. Let us consider the same parameters as in |l()j . namely 7 = 5/3, /3 = 1, the 
same abundance (^ = 0.619, ^c = 1.18, /-in — 1.43) as well as assume tsf — ^cooi and use l(TT|l for 
calculation of Atcooi- 

If shock wave moves in the uniform medium, then - with the use of Eq. I|18|l - the transition 
time is 

itr = 2.84 X 104£;^/"n-9/i7 yr (31) 

where E51 = i?sN/(10'^^ erg). The gas element which first cools (at tcooi) was shocked at ti which 
follows from Eq. flt^ : 



ti = 3.67 X lO^E^i 



4p4/17 9/17 



yr- 



(32) 




Figure 2: The evolution of the deceleration parameter m and different reference times for the shock 
motion in the uniform medium. Solid line - numerical calculations ^], thick dashed lines - Sedov 
solution (till th) and analytical solution ^^ (after Tgf). The dimensionless reference times are 
Tsag = 0.654, Tdyn = 0.802, Ttr = 0.855, n = 1.01, ri = 1.10, m = 1.57, now = 0.047, m = 3.03. 



The function m{T) reaches his maximum at radiative phase at t„ 



6.18 [ini. 



The time of the shell formation is given by Eq. ((SJ: 



isf = 5.20 X 10^Eti'\^'^" 



yr, 



(33) 



so that tsi/tt-c — 1.83. The time when the radiative losses of the shocked gas reach their minimum 
is ®: 

ilow = 1.60 X 103T3-J/^£;5\/\-1/3 yr (34) 

where T307 = Tiow/(3 x 10^ K). Under assumption that radiative losses does not change the shock 
dynamics till thj, with the use of Sedov solutions for the shock motion one have from Eq. Q that 



ihi = 1-04 X lO^Tj^ 



5^-5/6 pl/3 1/3 



o5 



51 



yr 



(35) 



where T2cb ~ rhi/(2 x 10^ K). The fluid temperature drops faster due to cooling than due to 
expansion from time 



idyn = 2.66xl0%\/^i-9/i^yr 



(36) 



The time when one may expect to have the temperature decrease downstream close to the shock 
is 



is 



2.17xl0Xi"n-9/"yr. 



(37) 



The Sedov solutions give at time itr the shock radius _Rtr = 19-E5( Uo pc, the shock velocity 
I?tr = 260 -Eg( nj km/s, the post-shock temperature Ttr = 0.95-10^ £'5( nj K and the swept 

up mass MtotiUv) = 10^ eH'^"^ Uo'^'^'^ Mq. 

The above reference times are shown on Fig. |2 together with evolution of the deceleration 
parameter m{T) calculated numerically ^D]- The analytical solutions for the adiabatic ^^1 and the 
radiative shock 1^ are also shown. Numerical result is found for supernova energy _Esn = 10^^ erg 
and interstellar hydrogen number density rio — 0.84 cm~'^. With these values, the times are 
isag = 2.4 x lO^yr, tdy„ = 2.9 x lO^yr, U, = 3.1 x lO^yr, t^ = 4.0 x lO^yr, i^f = 5.7 x lO^yr, 
iiow = 1-7 X 10'^ yr, ihi = 1-1 x 10^ yr; the intersection time is t\ = 3.6 x 10** yr [16j. The function 
raij) reaches his maximum during the radiative stage at imax = 2.3 x 10^ yr JS]. Results on Fig. [21 
are presented in terms of the dimensionless time t — t/t because the analytical solutions allow for 
scaling (numerical results for various input parameters differs by oscillation transient only; see e.g. 
Fig. 8 in |10p. The dimensional scale for time determined from fitting of analytical and numerical 
results \si = 3.6 x 10^ yr P,6 . 

It is apparent from Fig. |21 that the transition time ttr is a reasonable estimation for the end of 
the adiabatic stage while igf could be the time when one can start to use the radiative solutions 
P^ coming from the PDS model of McKee & Ostriker ^2- The duration of the intermediate 
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Figure 3: Numerical calculation of evolution of the deceleration parameter m from JO] (thin black 
line) and from [Hj (thick gray line). The transition and shell formation times from are marked 
by "C". 



transition subphase is {tsi — Ttr)/Ttr = 0.83 times the age of SNR at the end of the adiabatic stage, 
i.e. almost the same as duration of the adiabatic stage itself. This means that there is a strong 
need for a theoretical model which describe evolution of SNR in this subphase. 

For estimation of reference times, a number of authors i^ 1211 El 1221 keep a bit different 
approach from that used above, namely they use the approximation of the equilibrium cooling 
function with (3 = 1/2 and the Kahn solution for cooling time (|16(l . Let us compare the results 
of this approach with those obtained above. The evolution of the deceleration parameter in the 
refereed approach is presented in 1^1. There is also the same definition of the time of the shell 



4,73/14^-4/7 



formation tgf — tcooi- The estimation is ts[,c ~ 4.31 x lO^i^g^^^'no^' ' yr for their abundance and 
the cooling function (|30|l . For the parameters used in the numerical calculations -E51 — 0.931 
and Ho = 0.1 cm^'^ the time is isf,c = 1-58 x 10^ yr while with the use of our Eq. H33|) we obtain 
isf — 1.73 X 10^ yr. The both estimations arc close. Analytical solutions shows that, before tn- and 
after tsf, the evolution of dynamic parameters of the shock can be expressed in a dimensionless 
form, i.e. independently of E51 and rio. The behavior of the shock velocity depends however on 
these parameters during the transition subphase; the difference is in the frequency of oscillations 
(Fig. 8 in ^Ol)- Nevertheless, as one can see from this figure, the strong deceleration of the shock 
right after itr up to the first minimum is almost the same for different parameters, i.e. can also 
be scaled. We use this property in order to find the scale factor t for calculations being done in 
1^1. Namely, the fit of curve TO(r) from |6| to that of 10 (within the time interval from itr to the 
first minimum) gives ic = 1-05 x 10^ yr. The both calculations of the transition to the radiative 
stage agree rather well as it may be seen on Fig. O The dimensionless times for results in 1^1 are: 
the shell formation time Tsf^c — isf,c/^C = 1-51 and the transition time (as it is follows from (|38|l ') 
Ttr,c = 7sf,c/l-92 = 0.785. Fig. |31shows that the both approaches for localization of the limits of 
the transition subphase - with the use of the nonequilibrium-ionization cooling function (|29|l and 
the simple estimation for Aicooi ()ll|l [TH] or with the equilibrium cooling function (|30|l together 
with Kahn solution for Aicooi (|16|l |S1 - give almost the same estimations. 



3.2 Shock in a medium with a power-law density variation 

Let us now consider the shock motion in the ambient medium with the power-law density variation 
p°(i?) = AR~'^. With the use of ^^, (dJ, ©, © and the definition igf = ^cooi one can show 
that the duration of the transition subphase is given by 



tsf 



'•cool 
ttr 



l + ?7 



v 



')/(!+')) 



(38) 



The shell formation time is always larger than the transition time itr, provided by the fact that 
77 > 0. The ratio 



tt 



f] 



1/(1+')) 



(39) 




Figure 4: The ratios of times for /3 = 1 (thick lines) and /? = 1/2 (thin lines) as it is obtained from 
l|55|l and 123. 




Figure 5: The ratio ai/Rn and adyn/-Rtr for /3 = 1 (thick lines) and /3 = 1/2 (thin lines). 



is also always larger than unity. Note that these relations do not depend on abundance and 7. 
The ratios between all other times may be found from H28|l . 

The consequence of times is tdyn < ^tr < ti < tsf (Fig.0J in nonuniform medium with increasing 
density. The time ii may be smaller than tn- and tdyn for the decreasing density medium. The sag 
time isag < Ui for w > —6 only. 

Fig. 0] shows the two ratios (|38|l and (|39|l as a functions of ui for two values of (3. Namely, 
the ratios ti/ttr ~ 1.3 and tsf/ttr ~ 1-6 -^ 1.8 are almost the same for shock in the medium 
with increasing density {uj < 0). Therefore, in case of a uniform medium and a medium with 
increasing density, there is a need of introduction of transition subphase with duration more than 
a half of SNR age at the beginning of this subphase, tn-. The transition time tn- and therefore 
the transition subphase tgf ~ ^tr oc ttr are less for higher density and lower initial energy: itr ex; 

p(2+2/3+w)/(5 

in numerical calculations (Fig. 8 in 10 ). 



^^.^.H-r-;/*^-(7+2/3)/5 ^^^^.^ S^ 11 + 6(3 -uj{5 + 2/3). Such dependence on density is also visible 



3.2.1 Medium with decreasing density 

It seems that the formulae (|38|) and H39I) suggest for the case of decreasing density that the PDS 
radiative stage can even begin right after the end of adiabatic stage: tsf/^tr — > with uj — > 
3(1 + /3)/(2 + /3). Another result, already stated in 0, also follows: there will be no radiative shell 
formation for uj > 3(1 + f3)/{2 + (3). In order to understand the reasons of such behavior let us 
consider more details. 

What is the coordinate oi of the element which cools first? This element was shocked at 
ti = r7i/(n-';)t^j.. The Sedov radius at this time is i?(ti) = ai = r72/((5-a;)(i+r,))^^^^ ^Yms the 
coordinate oi > i?tr if w < 1.4 {(3 = 1) as it is shown on Fig. [3 The ratio ai/Rn is close to unity 
and is almost the same for such uj, i.e. the fluid we are interested in will be shocked soon after itr- 



However, li u; > 1.4 then ai — > quickly with increasing of a; from 1.4 to 2, i.e. the element which 
cools first is already inside the shock and may be in a very deep interior. The situation looks like 
that there could not be any "radiative shell" in a common sense. 

It is clear that the trend ist/^tr — > does not mean that radiative processes in the shock develop 
quickly for lo > 1.4. The transition and the shell formation times correspond to different processes: 
itr comes from comparison of the initial thermal energy density of the shocked fluid with radiative 
losses though tgi — tcooi is a time when the first cooled element appears. The two mentioned 
processes have place in vicinity of the shock if ambient medium is uniform or with increasing 
density. Numerical results suggest that they may be used for approximate estimates of the limits 
of the transition subphase in such media. However these two process are separated in space for 
media with decreasing density. It could be, that one (or both) of the times itr and igf may not be 
suitable to mark stages of SNR in medium with decreasing density. 

The cooling of shock moving in the medium with decreasing density differs from a commonly 
accepted scenario of the "thin dense shell" formation and should be studied in more details in the 
future. 

4 Conclusions 

The common approximate scenario of SNR evolution consists of the free expansion stage, the 
adiabatic phase and the PDS radiative era. It is shown that it is necessary to consider also 
additional subphase between the adiabatic and the radiative stages because this subphase lasts 
more than half of SNR age it has at the end of adiabatic stage. 

The analytical estimations on the ratios between the reference times which characterize the 
transition of adiabatic SNR into the PDS radiative stage - ttr, tgf and ti - does not depend on the 
initial parameters of SNR and IMS (energy of explosion, number density in the place of explosion, 
7 etc.) except of the density gradient (i.e. uj) and assumed /3 which causes rather small effect. This 
result is also visible in the numerical calculations for case of the uniform medium (Fig. 8 in |ll)|): 
except of the oscillations (which is indeed different for different rio) the durations of the transition 
subphase in terms of the transition time are almost the same for different values of ISM density. 

The ratio isf/itr ~ 1.6 for shock running in media with constant or increasing densities. The 
transition time however depends on the energy of explosion, the density of the medium and the 
density gradient: ttr oc E^j^'A~'''-^^ with a > and b > for shock in a medium with po oc i?~". 
This means that the transition subphase is longer for higher explosion energy and smaller density. 
The dependence of itr on this parameters are stronger for higher lo because the functions a{uj) and 
b{ui) increase with uj. 

The hydrodynamical properties of the shock in media with a; > seem to cause a trend to 
absence of the radiative phase in a common sense. The cooling of such shocks differs from a 
commonly accepted scenario of the "thin dense radiative shell" formation and should be studied in 
more details because it is important for models of nonspherical SNRs which could be only partially 
radiative. 
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Appendix 1. 

Approximation of the temperature evolution in a given fluid 

element downstream close to the strong adiabatic shock 

In order to simplify the estimation of tg^g and tdyn, let us approximate the distribution T{a) — 
T{a,t)/Ts{t) downstream close to the strong adiabatic shock; here a is Lagrangian coordinate, 
T = T/Ts and a = a/ R. Note, that hereafter in this Appendix we use the normalized parameters, 
i.e. divided on their values on the shock front; thus we skip the overlines in the notations. We are 
interested in the approximation in the form 

T(a) « a^''^'^'"^ (40) 



The value of k is given by 

/ d\nTia)\ , , 

where T(a) is the profile from Sedov |19| solutions. The equation of the mass conservation and 
the equation of the adiabaticity applied for the case of the shock motion in the medium with the 
power-law density distribution give the distribution of temperature T(a) — P{a)/ p{a) ,29, 

T{a) = (^)'"' «^^-'+" {riafrMV^' (42) 

where r is Eulier coordinate and ra = dr/da. Instead of Sedov profiles for r(a) - which is quite 
complex - we use the approximation 

r{a) = a(^-i)/T exp (a(a^ - 1)) (43) 

where a, (3 are constants; this approximation gives correct values of r and its derivatives in respect 
to a up to the second order on the shock [221 • Substitution H41|l with (|42|) , H43|) and with expressions 
for a, /3 from |29j yields 

_ 2(8-(7 + c.)(7+l)) 

(7 + 1)2 • ^^^^ 

For 7 = 5/3, k = 1 - 3w/4. 

The approximation (|4()|l underestimate Sedov temperature. The smaller a the larger difference. 
It is about 20% at a « 0.5 (that corresponds to r « 0.8). 

Appendix 2. 
List of times 

isag "sag" time ^, radiative cooling begins to affect the temperature distribution downstream 

of the shock; 
idyn "dynamics-affected" time T, the temperature of a fluid element shocked after this time 

decreases faster due to radiation than due to expansion; 
ttr "transition" time ^01 , estimation of the time when the deviations from Sedov solutions are 

prominent; Sedov solution may be approximately used till this time; 
Atcooi "cooling" time |23[I2U[0T] . a shocked fluid element cools during this time; 
is "shock" time [2011^1 moment when the shock encountered given fluid element; 

ii moment when the shock encountered the fluid element which cools first |20II21| : 

ic sum of is + Aicooi; 

icooi "SNR cooling" time |20[ 1211 1^ . the minimmn of ic, i-c the age of SNR when the first 

cooled element appearsj 
isf "shell- formation" time |20ll21ll^ . approximately after this time the shock may be described 

by the radiative PDS model; 
iiow moment during the adiabatic stage when the radiative losses of the decelerating shock wave 

reach their minimum value; 
ihi moment when the radiative losses of the decelerating shock wave reach their maximum value 

ipiini; 

ii "intersection" time 16 , moment when two functions - adiabatic R{t) and radiative i?sh(i) 

intersect; 
imax moment during the radiative stage when the function m{T) reaches its maximum |lti| : 
i timescale; 

r dimensionless time, r — t/t. 
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